function q = ElnqZ(vars)

% E[ln q(Z)]

Z = vars.Z;

q = sum(sum(Z.*log(Z+realmin)));
